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ABSTRACT 


This thesis presents an analytical technique for the 
design of compensation for multi-loop servomechanisms. A 
completely general multi-loop system is represented ina 
standard block diagram. For this system a standard deter- 
minant is formed which uniquely represents each block, loop 
and path in the block diagram. Based on this determinant, 
the effect of all possible compensators on the coefficients 
of the system characteristic equation is found. This infor- 
mation is presented in a Location Matrix. A model character- 
istic equation is synthesized from the static and dynamic 
performance specifications. Comparison of the uncompensated 
characteristic equation to the proper model equation shows 
which coefficients must be altered. The Location Matrix 


will show all the compensator-path combinations which will 
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affect the desired change. The method requires no intuition, 
and so can be applied by persons without extensive training 
in the servomechanism field. A step-by-step procedure for 


application of this technique is included as a summary. 
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OBJECT 


Ihe object of this thesis is to present an analytical 
technique for the design of multi-loop servomechanism 
compensation. This technique is set forth under the 
following topics: 
(a) Formation of the Characteristic Fquation. 
(b) The Effect of Performance Specifications 
on the Coefficients of the Characteristic 
Equation. 
(c) The Effect of Compensation on the Coef- 
ficients of the Characteristic Equation. 
(d) The Choice of Compensation. 





Chapter 1 


INTRODUCTION 
and 
METHOD OF APPROACH 


l.l Introduction. The problem of compensation design for a 
multi-loop servomechanism is not a new one. Indeed, most 
standard servomechanisms or control systems texts give it 
some mention. Most such works treat the problem as an 
extension of the ee problem, and attempt to apply 
single-loop concepts to its solution. 

The penalties of such an approach are numerous and 
serious. The first such penalty arises from the technique 
of getting the system into a single loop configuration. This 
involves a block diagram reduction or some similar procedure. 
As a result of this reduction, the identities of the original 
blocks, and the physical components they represent, are lost. 
Far more serious than this, in a compensation design problem, 
is the fact that the information about the arrangement of the 
components is lost. The analysis of this reduced system would 
show algebraically what kind of compensation was needed. It 
would be impossible to transform this information into a given 
compensator fed around certain main path elements. The whole 
sense of element and path configuration no longer has any 
correlation with the algebra. | 

There are two ways out of this dilemma. The designer 


can insert all possible compensators, in all usable paths, 
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and re-do the analysis until he finds one combination that 
gives the desired result. Obviously this is a monumental 
task for any system of reasonable complexity. There are 
graphical techniques which eliminate some of this work, but 
they still involve trial and error. The alternative is the 
intuitive approach. Here the designer must have such nen 
experience that he generally knows the effect of all com- 
pensators, and so can eliminate from consideration many 
which could not work. 

In this thesis an analytical technique is presentec which 
requires no intuition. It examines the effect of all possible 


compensators without the accompanying drudgery. 
1.2 Method of Approac 


A general multi-loop control system, consisting of all 
possible direct, feedback and feed forward paths, is formulated. 
Its components are presumed to be describable by ordinary, 
linear integro-differential equations, with constant coef- 
ficients. When inputs to this system are mentioned, the 
terms "step", "ramp", and "parabola" of input will be used. 

A step of the input quantity is that input which yields zero 
error when there is unity feedback from output to input. For 
a Type I positional servo, a step would be a step of position. 
For a Type II velocity servo, the step would be a step of 
velocity, etc. By treating such a general system, the ensuing 
results can be applied to nearly all practical systems. 

This analytical compensation design technique uses as its 
cornerstone the characteristic equation of the system. After 


a brief discussion of the methods for choosing direct path 


components, a method for forming this characteristic 
equation is described. The system is represented ina 
standard block diagram. Using circuit analysis techniques, 
the interaction of the blocks, path and loops in this dia- 
gram is shown in a characteristic determinant of the system. 
It is from this determinant that the characteristic equation 
is written. 

The next step in the procedure relates performance speci- 
fications to the coefficients of the characteristic equation. 
Static specifications such as "following error" and output 
in response to a load torque, place limits on coefficients 
or on ratios of coefficients in the characteristic equation. 
Because of the vast amount of work which has been done on 
simple second order systems, most dynamic specifications 
are expressible in terms of the parameters of such systems. 
The compensated system is made to have a response not sig- 
nificantly different from the specified second order response. 
This is done by making the quadratic characteristic equation 
of the second order system a factor of a model characteristic 
equation for the compensated system. All other roots of this 
model equation are forced to be so large that the second order 
roots dominate the response. 

The effect of all possible compensators on the coeffi- 
cients of the characteristic, equation is found through the 
use of the characteristic determinant. All possible paths in 
the system are uniquely represented in this determinant. 
Closing any of these paths with a compensator adds a cofactor 


to the existing determinant. A rapid method for determining 





the effect of this path closing on the coefficients of the 
characteristic equation is presented. This information is 
displayed in a Location Matrix. 

At this stage of the design process, the designer has a 
model characteristic equation representing the specifications, 
and an actual characteristic equation representing the uncom- 
pensated system. A comparison of these shows which coefficients 
must be altered by compensation. The Location Matrix indicates 
what compensator, in what location, will effect the required 
change. 

This thesis confines itself to the problem of assuring 
dominance of a desired second order response. The technique 
itself is not so limited, and could be used to force a system 


to perform like any desired third or higher order systen. 





CHAPTER 2 


BASIC DESIGN CONSIDERATIONS 


iv Tnarroduetion 


Some servomechanisms perform a computing function. 
Examples of this type of servo are found in fire control 
systems. In these servos, the output or outputs may be com- 
plicated mathematical functions of several inputs. The mathe- 
matical relationships to be instrumented will dictate the 
choice of equipment to be used. The type of servo of primary 
interest in this thesis, however, is one which drives a load 
so that its action duplicates the action represented by the 
input signal. An example of this. type of servo is the drive 
system for a gun. Its input is a voltage or shaft rotation 
which represents the desired orientation of the gun barrel. 
It is the function of the servo to control the power that 
drives the gun mount so that the desired gun barrel orien- 
tation is achieved. The preliminary choice of equipment to 
perform this type of function is considered in this chapter. 
These A result in a basic system with sufficient 
power and accuracy to meet the fundamental requirements set 
by the specifications, Subsequent chapters deal with modi- 
fying the basic design to meet specifications concerning the 


more detailed nature of the servo's response. 





ex Components of & Typical Servomechanism 


Figure II-1 is a diagram of a typical servomechanism. Its 
purpose might be to control the angular position and/or velocity 
of a shaft. The shaft is represented by the block labeled 
"Load". The input oT command might be the angular position 
of a shaft from a computer. The command receiver and the out- 
put measuring device convert the input and output into quanti- 
ties which can be easily compared. The comparator subtracts 
these quantities and produces a signal proportional to the dis- 
parity between input and output; 1.e., error. Since this error 
may be a small quantity, it must be amplified. It is then’ 


applied to a controller which actuates the prime mover. The 


prime mover drives the load through the gear train so as to 


null the error. 





Figure II-1 


A servomechanism can be considered to consist of paths 
through which information flows. Those paths which conduct 
information in the direction from input to output are termed 
forward paths. Those which conduct in the opposite direction 


are termed reverse paths. In essence, the forward paths of a 





servo transform the input into the output. The reverse, or 
feedback, paths insure the proper degree of correspondence 
between input and output. 

In a relatively simple system, there is usually one for- 
ward path. However, other forward paths may result from com- 
pensation considerations. To avoid ambiguity, the forward 
path which performs the primary function of changing input to 
output will be designated the main path. Other forward paths 


will be designated feed-forward paths. 


3. Some Considerations Governing Choice of Components 


The need for a servo implies that the element which generates 
the input has insufficient power to control the load. Since 
the purpose of the servo is to supply the needed power, choice 
of the prime mover is a logical starting point in selecting 
main path components. The characteristics of the load will 
be known. These characteristics may be expressed in terms of 
its inertia, the viscous friction to which it will be subjected, 
external torques which may act on it and the like. The speci- 
fications will also state the velocities and accelerations re- 
quired of the load. This information will set a lower limit 
on the power and torque capacity required of the prime mover. 
The power supply available will further limit the choice. Other 
requirements on the prime mover will arise from cost, weight 
and size limitations, life expectancy, etc. With these limits 
in mind, the designer must then select a prime mover from those 
commercially available or else design his own. 


Choice of the prime mover will normally place limits within 





which a gear ratio must be chosen. Ahrendt! discusses the 
war ious) factors which influence this choice. Such factors 
as low cost or smoothness of output at low velocities may 
place conflicting demands on the choice of gear ratio. Also, 
the inpedance matching aspect of gear ratio selection may be 
of great importance. In order to choose the optimum gear ratio, 
the designer must have a thorough-going knowledge of the speci- 
fications as well as of the characteristics of the prime mover. 
Choice of the prime mover will probably narrowly limit 
Chewec of behe controller, If a variety of controllers is 
available, the considerations of fineness of control versus 
cost, weight, etc. must be brought to bear. 
Insufficient information will be available in the speci- 
fications to allow choice of the amplifier to be used. The 
gain to be provided will be determined by compensation considera- 
tions. This will not be known until later in the design process, 
When it has been determined, the amplifier may be chosen, In 
addition to the conditions imposed on amplifier selection by 
such factors as available power supply, noise generation, and 
operating environment, the designer must bear in mind the impor- 
tance of linearity. Although sufficient band width is not diffi- 
cult to obtain, saturation effects may present a problem, It 
is uneconomical for a system's amplifier to saturate before its 
prime mover does, as this means unusable power has been purchased. 
Ihe choice of comparator, measuring device and command 
receiver will be a question of obtaining the best linearity and 
accuracy within the limits set by cost, weight and size speci- 


fications. These devices are described in detail in several 
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servomechanism tome 


l. Determination of Component Transfer Functions 


Having chosen the main path components, their transfer 
functions must then be determined. For components which have 
their output independent of the frequency of the input, this 
merely involves determining their sensitivity. It is desirable 
that this sensitivity will be constant over the range of interest 
in the problem, If sensitivity is not constant, it must be 
assumed so in the remainder of the problem by assigning it a 
representative value. 

The transfer functions of frequency dependent components 
may be determined mathematically by application of basic a-c 
circuit theory and the Laplace Transformation. This technique 
is described in detail by Tee If the electrical parameters 
of the component are unknown, it may be necessary to conduct 
tests on the equipment to determine its transfer function. 

These may take the form of frequency or transient response tests. 
The results of these tests will permit derivation of the trans- 
fer function, Examples of this method are presented in refer- 
ences 5 and 6. 

Having obtained the transfer functions of the equipment 
to be used, it is possible to express the system in this stage 
of its design by a signal flow chart or standard block diagram. 


EU SNpHUSENOÍ design is disieussed in the next chapter. 
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CHAPTER 3 


FORMATION OF THE CHARACTERISTIC EQUATION 


3-1 Introduction. The characteristic equation is an equation 
in s domain which represents the integro-differential equa- 
tions in time space which describe the flow of information 


through a system.’ 


The roots of this equation allow the 
determination of the transient performance and the stability 
of the physical system so described. The purpose of this 
chapter is to develop a standard method of forming this 


characteristic equation for a multiple-loop servo system. 


3.2 Standard Block Diagram. A set of transfer functions are 
developed, as described in Chapter 2, for elements of the 
servomechanism. The block diagram portrays the interconnec- 
tion of these various elements to form the overall system. 

The very nature of the block diagram assumes that the transfer 
function of the tandem conbination of two blocks is the product 
of the individual transfer TUS." The overall perfor- 
mance of the system is found by combining these blocks 
according to the standard rules for their manipulation. 
Unfortunately for the designer, the combining of the blocks 
destroys their individuality and the identity of the mechanism 
or circuit they represent. What is needed is an alternate 
representation of the system which pictures it in more detail 


than the reduced (combined) block diagram. This could take 





ie 


a new form, as in the signal-flow diagrams originated by 


1 or it could be a more ordered manipulation of 


10 


o. d. Mason, 
the block diagram, as proposed by Yaohan Chu. It is to 
the latter method that we now devote our attention. 

Certain advantages accrue from standardizing the block 
diagram. In this regard the following symbols are defined: 

a. The Nodal Point or Node: These are summation points 
in the main path of the system. Fach node may receive a 
number of signals, but delivers only one. The nodes will be 
lettered in sequence, a, b, c, . . n, from left (input) to 
right (output). 

b. Ihe Pick-off Point: This is a point at which there 
is one incoming signal and two or more outgoing signals. The 
outgoing signals are identical with the incoming one. 

c. The Block Transfer Function: The transfer functions 
in the blocks will be designated by double subscript notation. 
This notation is defined for the three types of blocks as 
follows: 

l) Direct Path Block: These blocks are noted Gy where 
the "n" refers to the node preceding the block. "N" is the 
upper case form of the same letter. 

2) Feedback Block: The feedback blocks are subscripted 
Gym‘ The "N" is the same as the upper case subscript of the 
direct path block from which they receive their input. The 
"m" is the designation of the node to which they feed. 

3) Feed-forward Block: Such blocks contain the symbol 


Gun‘ The lower case subscript refers to the main path node 
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preceding the block. The upper case subscript is the same as 
the upper case subscript on the main path block which feeds 
into the same node. 

As examples of the above, consider the following three node 


system: 


Figure 3-1 


cC Output 





(3; (9, and (©) are nodes. u symbolizes a pickoff point. 


Although the system above is represented as having but one 

input and one output, there is no limit to the number of 

either. Internal noise generated in one block may be considered 
as an external input to the node following that block. If 

there is more than one input, each is designated by a Roman 
numeral. If there is a transfer function between an input and 


its entering node it will be subscripted as follows: 


(Input "T" G 


m e 

E um | 

— MA 

o o ^ 

E g - E 
[3 f E 
"- A 
= @ 
Ë = Lb 
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can obtain a set of simultaneous equations for the system 
represented by a standard block diagram s applying the network 
analysis technique of Guillemin.2? This entails writing the 
equation for the signal summation at each node. If we iden- 
tify the signal out of a node with the lower case letter of 
the node, we have the following equations in our example: 

At Node a: 


a = 1I(G,,) - 0(G ¿Gp,) - a(G,¿Gp,) - DGypGg,) 


At Node b: 
b= a(G, 4) - b(G, Gap) - c (Go cEcp? 


And, at Node c: 


c 7 a(G,g) * b(Gyy) - c(G 6G) 


Rearranging these equations we have: 


(1 + Carpa)? K (OG, gg b h (Goclca?e 7 (Gia) I 


"d 


)b + (G 0 


re FG se 


bBUBL 
(-G,p?a + (-G,5)b + (1 + Gack = © 


Inspection of these equations reveals a symmetry of arrange- 
ment which would permit writing the equations almost without 
reference to the block diagram. All the forcing functions are 
on the right. On the left are terms which contain components 
and signals within the system. We can express these equations 


in matrix notation as follows: 


LI HA E = 
j 
LJ 
A A 
- 
- "C 
u ' 
=) 
Es 





Tp 


LR UU Cop Ba CGocóca a 1(6,,) 
Goa 1 + G Bfsp Cecêce PALOS 0 
“Gap “Go l * Cocfco]| e 0 


We call the determinant of the matrix of the coefficients 
the characteristic determinant of the system. To this 


determinant we assign the symbol A- 


3.3.2 Manipulations of the Characteristic Determinant. 
Independent of its physical significance, the determinant of 
the coefficients of a set of linear equations has certain 
properties. These properties are illustrated in the tech- 
nique of equation solving with determinants and in the rules 
for expansion of a determinant. 

Ihe solution of the equations is accomplished as follows: 

1) Replace the column corresponding to the desired 
unknown with the column of forcing functions. The resulting 
determinant, divided by /\, is the value of that unknown in 
response to these inputs. 

2) Since the determinant represents a linear system, the 
principle of superposition is valid. Hence, in the case of 
multiple inputs, the output resulting from each may be com- 
puted separately and the individual results added. 

3) The output obtained by replacing the last column of 
the determinant with the input column is the signal from the 
last node. It must be multiplied by the performance function 


following the last node to obtain the system output. 





16 


The characteristic (homogeneous) equation describing the 
system is obtained by setting the characteristic determinant 
equal to zero. 

The normal rules for expanding a determinant apply to the 
characteristic determinant. In general they require the expan- 
sion by co-factors of some row or column. Since we will 
usually be interested in the final output, we should expand 
by elements of the last column. To the co-factor of the ith 
row and jth column we assign the symbol Ars Similarly the 
element in the ith row and jth column is called G4 4* Using 
this notation, the expansion of an "n by n" characteristic 


determinant by elements of the last (nth) column is: 


n 


=, Cin Mn 


A 


As an example of the procedure described above for finding 
system output in response to a given input consider the 


following: 


Input £o node "a" is I, Gra = 1 


=G I 726.26 


aA bB Bb 


(Goo) I sa g 
aB bB 


SS å 


Output 


Output * I(G, C) [G, Gy * Gap * GS rpg 


3.3.3 Properties of the Characteristic Determinant. In 


addition to properties inherent in all determinants the charac- 


teristic determinant has features peculiar to the system it 





n 


represents. 

1) The characteristic equation, obtained by equating /N 
to zero, is the sum of unity plus a number of closed loop 
performance functions. By the nature of the elements which 
comprise the determinant, these functions will be ratios of 
polynomials in s. When this sum is placed over a common 
denominator, its numerator will oe the usable form of the 
characteristic equation of the system. 

2) Each element of the characteristic determinant 
corresponds to a path in the standard block diagram. 

3) The diagonal elements of the characteristic deter- 
minant correspond to the paths of inner loops of the nodes. 
The elements above these diagonal elements are all feedback 
paths. Those below the diagonal are feed-forward or direct 
paths. 

4) Each feedback and feed-forward performance function 
appears in only one element of the determinant. Each direct 
path performance function may appear in more than one element. 

We now have at our disposal a standard block diagram. 
Ihis block diagram maintains the integrity of each element 
comprising the system it represents. We form this block 
diagram to aid us in compensating the system. To this end 
we need only include those nodes or pick-offs which represent 
points in the system physically usable for compensation. All 
performance functions between such points should be combined 
into one block to simplify the ensuing algebra. The expansion 
of the standard determinant representing the block diagram gives 
the characteristic equation of the system. Chapter 4 will relate 


performance specifications to this characteristic equation. 
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CHAPTER 


—— SS ee ee eS nO 


COEFFICIENTS OF THE CHARACTERISTIC EQUATION 


le Introduction 

The performance specifications state the desired system 
output in response to actuating or disturbing inputs. The 
output is specified in terms of its steady state value and its 
transient behavior. | 

The characteristic equation which was evolved in the pre- 
vious chapter describes the system performance in the s domain. 
The standard analysis approach is to find the roots of this 
characteristic equation and transform it to the time domain 
where actual and specified performance can be compared. 

The purpose of this chapter is to relate the specifications 
directly to the coefficients of the characteristic equation. 
This eliminates the need for finding the roots, plotting the 
response, and other laborious analysis techniques. 

2. Static Performance Specifications 

The static performance specifications describe the steady 
state value of the system output. In response to an actuating 
input the specification will normally state an allowable error. 
In response to a step we would expect no error. A following 
error will be stated for the response to a ramp input. A 
maximum allowable output is usually stated for the response to 
a disturbing input. These disturbances might be noise or some 


form of load torque. 
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Evaluation of the steady state output to a step, ramp or 
parabola of input position is most useful in studying the effect 
of disturbance signals. As established in Chu's Ti tele omn 
the performance function relating the input at any node with 
the output at any other node is 

Output _ yA) (1.3) 
nput PX 
where: iO node into which input is fed, 
jo node just preceding output block, 
Gg a performance function of output block, 


th 


As; © cofactor of the i” row and the jee column 


of the standard determinant. 


The output in response to a given input, I(s) is 


——— (4.2) 


The steady state value of this output is obtained by applying 


the final value theorem as follows: 


I(s) G 
Output s = Limit ^? ET (4.3) 
=> 0° 
E anis /\ = Characteristic Equation (loli) 


r 
a 
where Dí common denominator of the ratios of 


polynomials making up the expansion of /\. 


s I(s) 6,,/\,;; DA 
Characteristic Eqn. (4.5) 


Output E m Limit 
3 s -— 0 
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For a step of disturbance, i.e. I(s) = A/s 


Output, = Limit E DEEP | since a 


s > 


O 


is the limit of the characteristic equation as s goes to zero. 


For a specified maximum output over input ratio, k, where 


$a output ss 
k > = le A; DA| (4.6) 
a 
O 


A close examination of the equation above will reveal the 
effects of such a specification on the a, coefficient of the 
characteristic equation. 


Case I - If k is a number not equal to zero, then a, is 


Specified as being greater than the Limit Ea PA 
k 


s = 0 
This may or may not limit the size of a, since this limit may 
be identically zero. 
Case II - If k is specified as zero then the 


Limit G CD must be zero since a, cannot be 
amio [es A, ^| o 
infinite. 

In Case I where the Limit |G, ; DA | is a constant 

S70 DM 
we have set one condition on a, which may be expressed in 
equation form as: 
ao ZN 


Note that the Limit ce. Ar ED may be constant containing 
S 0 JJ c 1) A 


variable gains of the uncompensated.system. The condition on 
a, may be met by setting these gains which determine the stiff- 


ness of the system. This will then place a lower limit on 





ex 


these gains and restrict their use in future compensation. 


In Case II, if the Limit |G CD 0 then an 
crees [sts BAJS o e 


s must be introduced into the numerator of A ; which does not 


cancel an s in its denominator; or an s must be placed in EDAM 


which is not also placed in the denominator of Ay This 
will raise the order of the characteristic equation. The former 
may be accomplished by introducing derivative feedback from 
any block before the output block, to the actuating input's 
node. The introduction of the excess sin Di is brought 
about by putting integration feedback from other than the out- 
put block to the input node to which the disturbing inputs 
are TN 

To illustrate the reasoning behind the statements above 
consider a four node system represented by the following stan- 


dard block diagram: 


me) 
ipput qd g 


(disturbance) 
— Output ) 
? 


ote TANE ae Pros 







dD Da 
5 “Goa al O O 
- d O “Gop 1 O 
O O -Goc 1 


For a disturbing input at node "c" we wish the output to 


e. 3 3 0. 
be zero; but have postulated that pemi [Sap As, P| E 
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We need to introduce an s into one of the terms of this express- 
ion which will drive the response to zero without affecting the 
response to the actuating input. 

The steady state value of the response to an actuating 


step input, A/s, is (by equation l.5): 


Limit A ean Da] =A 


s>0 a 
O 


Since Gap and N are common to both expressions a change in 
either will affect both responses. If the s is introduced in 
PA it must also appear in the denominator of Sl so that 
its effect is not felt in the case of the actuating input. It 
must not, however, appear in the denominator of Ash: Thee2nd 
and lth rows of NM are contained in both cofactors so they 
cannot be used. The 3rd row appears in Ay) but not in PN" 
The only locations usable for an integration feedback which 
would have the desired result are positions in the 3rd row. 
In our example the only such feedback position is around Goc’ 

The other method of introducing an s in the expression for 
steady state output in response to a disturbance, without alter- 
ing the desired actuated response, is to put an s in the numer- 
ator of Az which does not appear in the numerator of Ary: 
This is accomplished by inserting derivative feedback in a 
first row position. The lst row is the only one in A3), which 
is not also in Air: 

In higher order systems there may be peculiarities of the 
determinant which would allow other positions to be used, or 
which would eliminate some of the above positions. Those 


explained above are the only positions possible if all terms 
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in the determinant are present. 

Computation of the output quantity is of prime import- 
ance in devising methods of making a system insensitive to 
disturbing inputs. In response to an actuating input we 
might better express steady state performance in terms of 
error since this permits an immediate and direct check on 
the accuracy of the system. If error is defined as the 
difference between the actuating input and its desired out- 


put, then the error equation is: 


E = Input - Output (4.7) 
This may be rewritten as: 
A . Output 
E h ei Input. (4.8) 


Expressing Output/Input as per equation(4.1) we have: 


l -G 
= I(s) (4.9) 


By the final value theorem, and again expressing A as in 
equation (4.4), we have 
Error... = Limit siI(s) Char.Eqn. - G. 5A; ; CD ^ (4.10) 


ss 
"ecc Char. Eqn. 


In response to I(s) = A/s: 


Error. = Limit A |Char Eqn - G,-/\,. CD (4.11) 
ss are Wei] AA 


Char. Eqn. 


Since the steady state error in response to a step will be 


defined as zero we can write: 
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Limit A [Char «Pan. - 83743 DA | 


> we 90 
pe Limit Char. Fan. (+.12) 
se 0 


Since Limit |Characteristic Equation] E, X 0, then: 


s> 0 


a, = Limit [eam CD A] (4.13) 


s> 0 
This defines another condition on ao? but one which will 
automatically be met if there is unity feedback from the 
output block to the actuating input node. 


In response to an input ramp, i.e. I(s) - Bye. 


Error- = Limit B Char. Eqn. - AA CDA | 


SS 
= Char. Eqn. 
(4.14) 
For zero error in response to a step, the Limit of the 
S0 
Characteristic Equation equals the Limit G CD =a. 
q q E 322843 ¡ÓN O 


Hence the constant parts of the numerator of equation 4.1) 


are equal. This means that s is a factor of the numerator. 


Error... = Limit B E tag. - 6.744; | (4.15) 
SE at 


S 


Error.. e 
B 


PME 
E 


E U ? 


O 


This equation puts a limit on the ratio of a,/a, in the 
ed characteristic equation. 

Tees ratio may not be compatible with the a/a, ratios 
of any of the model equations formed from dynamic specifi- 


cations. The restriction on a7/a, mày be relieved by makifig 
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the ramp following error identically zero, i.e., raise 
System to Type II. 

All the conditions on equation order, compensator loca- 
tion and coefficient magnitude which arise from consideration 
of static specifications should be held for simultaneous 
solution with similar conditions yielded by the dynamics of 


the problem. 


C. Dynamic Performance Specifications 


The dynamic performance specifications describe the 
desired transient performance of the system. Since they are 
meaningless in an unstable system, they all imply stability 
even if they do not state it. 

As opposed to the static specifications, where there 
is some uniformity, there are nearly as many ways of express- 
ing a desired transient behavior as there are Process Engin- 
eers. The following summary ^, while not intended to be all 
inclusive, lists the most common specifications and discusses 


how to calculate them. 


THE TEN MOST COMMON DYNAMIC PERFORMANCE SPECIFICATIONS 





NAME TYPE METHOD OF COMPUTATION 
1 
Gain Stability specifi- Compute maximum gain 
Margin cation 


for stability by Routh- 
Hurwitz criterion. Ratio 
of maximum stable gain 
to actual gain is gain 


margin. 





Phase 


Margin 


M Peak 


Frequency domain 
stability specifi- 


cation. 


Frequency domain 
stability spec- 
ification. 


26 


On Nyquist diagram, the 
angle between curve and 
negative real axis at the 
unit circle. On Bode 
plot, find a, (frequency 
of unity gain), compute 
phase shift at Was and 
subtract from 180 deg. 
For a well behaved system, 
a phase margin of 5 deg. 
means one overshoot in 


response to a step. 


The M curve is the mag- 
nitude of the closed- 
loop response vs. 
frequency. It can be 
calculated from open-loop 
transfer function or from 
values picked off Nyquist 
plot. M peak is defined 
as the peak of closed- 
loop frequency response 
normalized to low fre- 
quency value. 


M = A 
LA 


Circles of constant M 


can be drawn on Nyquist 





4 Stability spec- 
Damping ification. 
Ratio, Y 

5 Stability 
Damping Specification 


Factor, Cw 
n 
or Decrement 


Factor 
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plot. The locus of the 
system must remain out- 
side given M peak circle. 
(An M peak of 1.4 usually 


means one overshoot) 


Can be calculated from 
closed-loop transfer 
function or found directly 
from root-locus plot. 
Damping ratio is defined 
as E in a quadratic term 
of the form: 

se p celos + ay * 

On the root-locus dia- 
gram, © equals cos 9 
where © is the angle 
between the horizontal 
and the vector directed 
toward the gain point 


on the locus. 


May be found directly 
from the root-locus. 
Damping factor is defined 
by the factors of a quad- 
ratic system, as was 
damping ratio; it is the 


real portion of the roots. 





6 
Percent 


Overshoot 


7 
Bandwidth 


Direct evaluation 
of system 


stability 


Frequency domain 
Speed-of -response 


specification 
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Where the imaginary por- 
tion of the root deter- 
mines the frequency of 
the damped sinusoidal 
response of a system, 
the real portion deter- 


mines the rate of decay. 


May be calculated from 
response of system to 
step function input. 
Defined as peak of the 
response to a step-func- 
tion input as a percent- 


age of final value. 


Usually defined as the 
frequency at which the 
closed-loop frequency 
response falls to 1/2 
its low (zero) frequency 
value. Bandwidth is 
available directly from 
the M curve of the system. 
An analogous definition 
would be the crossover 


frequency of the system. 





8 Speed-of-response 


Rise Time specification 


9 Overall performance 
"Settle-Out" Specification 
Time or 

"Synchronizing" 


Time 
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One of the simplest 

definitions of rise time 
is l/bandwidth, another 
is 1/oy b odis Rise time 
is also defined as time 


to first zero of error. 


Settle-out time can be 
calculated analytically 
by the inverse Laplace 
transform, or by picking 
values off the root locus 
plot. This specification 
is most convenient if an 
analog computer is avail- 
able. Settle-out time 

is defined as the time 
required for system 
response to a step-func- 
tion input to approach 
and remain within a 

given tolerance of the 
final value. It is 
usually necessary to 
specify maximum over- 
shoot and steady-state 
error in addition to 


settle-out time. 





10 

ITAE or Integrated Overall 

Value of the Performance 
Product of Time Specification 
and Absolute 


Value 
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The ITAE specification 
is defined by the integral 


00 
IT" Sela 
0 


Minimize the value of 
this integral for optimum 
system performance. This 
specification yields a 
number that depends on 
system parameters. It 
would probably be an 
improvement to generalize 
it by normalizing to the 
rated output quantity. 
The dimensions would then 
be see”. Another poss- 
ible generalization 

would be to integrate 

to some arbitrary large 
number, such as 10 time 
constants, rather than 

to infinity. This would 
permit the specification 
to be applied to responses for 
which the system yields 

a finite final cretus 


ITAE is one of several 
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attempts to define an 
overall figure of merit 
for system operation. 
Rather than simply sum- 
ming the total error of 

a system, the error is 
progressively weighted 
more heavily as time goes 
on. This puts a premium 
on rapid, accurate settle- 
out, and allows for 
unavoidable initial 


large error. 


Due to the complexity and variety of responses of sys- 
tems of order higher than second, very little information 
is available in the literature concerning the characteristics 
of their response. However, the response of a second order 
system is easily computed, fairly simple, and can be Charac- 
terized by two parameters. If the system has an oscillatory 
response, the parameters most commonly used are € and On) 
the damping ratio and undamped natural frequency respective- 
ly. The use of these two parameters is quite natural since 
the coefficients of the system quadratic equation are simple 
functions of X and wo 

As seen above, dynamic specifications will not always 


be given in terms of č and E but translation to this form 


— p 
— Á— d 
— = 
Lr m 
m 
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is possible. As an example, Instrument Engineering, Vol. 1113 
provides graphs of calculated second order response in both 
the time and frequency domains which may be used for this 


14 


purpose. Chu gives several rules for translating time 


domain specifications into & and o and Thaler gives others.” 
Thaler also includes simplified graphs giving relationships 
between frequency domain specifications and X and wp as 
well as rules of thumb relating phase and gain margin to 
height of resonant peak. 

Translation of specifications into © and w OT relatable 
parameters may be other than precise. Analytical relation- 
ships may not exist, as in the case of phase and gain margins. 
Specifications may put varying đemands on the physical com- 
ponents of the system such that time dependent coefficients 
are required. They may prescribe a performance unrealizable 
by an order lower than third. In cases such as these, 
reversion to trial-and-error design techniques may be dic- 
tated. If the response of a higher order system can be made 
to closely approximate a second order system, then adequate 
correlation between response and the €t, o type specifications 
is possible. It is to the problem of assuring dominance of 
a given quadratic factor that we now turn our attention. 

A quadratic factor is said to have "dominance" if the 
response to the higher order system closely approximates that 
of the specified quadratic. A high "degree of dann e 
means a close approximation, i.e. small deviation of actual 


from second order response. The use for which the system 


is intended and the strictness of the specifications will 
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indicate the degree of dominance required. 
In an oscillatory second order system, the correlation 
between coefficients of the characteristic equation and 


the X and oO, specifications is 


ge + 20, S + w =00 (4.16) 


If our system has a cubic characteristic equation with 
one pair of oscillatory roots it will appear as follows on 


the s-plane: 





Figure 4-1 


where 23» the dominance factor, is defined as shown in 
figure 4-1. The model characteristic equation in terms of 


ai» © and w is: 


2 a. | 
(sta, )(s +2to, st, )=0 


(2% ro, Ja S r(2La,+1)a ^s*a aj? 20 (4.97) 


The degree to which the response of this system approximates 
a second order response depends on the size of 04 * A detailed 


examination of the effect of a, on the deviation of a third 
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order response from that of a second order is shown in 
Appendix A. 


A system with a fourth order characteristic equation 


might take either of the following s-plane configurations: 





Figure +-2a Figure +-2b 


The model characteristic equation in the case shown in 


figure 4-2a is: 
(sta, )(sta,w_ )(s+2Lw stw “) = 0 
i m 2n n n 


; s*« (2c*9, Ja, sJ« «2C p 8.) ^s? (9, *2t8, Ja 3sep, n," = 0 


(4.18) 


For the case in figure hà-2b the model characteristic 


equation is: 


2 Econ tS E 
(s +20, € 70, Sto, o )(s +20, w Sto, )=0 


s* (2€, *2a, C), s) (19a, C Cota, Ja 252 (2a, C, *2t a, a, s 


2m. s 0 TIE 


ta 
qe 
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In choosing the a's for the case in figure 4-2a a 
reference to the curves of Appendix A is helpful but not 
exact. Each a must be sufficiently large so that its 
effect alone would be within the tolerances allowed. A 
conservative choice of the individual a's is necessary to 
insure against the possibility of accumulating error of the 
same sign. Determination of the effect of more than one 
real root is a fruitful area for further work. 

In the case where there are two pairs of oscillatory 
roots, there is a slightly different problem. Here it is 
not only the decay of the non-dominant pair, but the size 
of the coefficient of their sinusoidal contribution to the 
ida response, which determines their effect. In the 
extreme case, where they are on the imaginary axis, this 
coefficient is approximately as times that of the domi- 
nant pair. -— 

The extreme condition described above also serves to 
place a lower limit on the size of the coefficients of the 
characteristic equation in order to have stability. If 
E = O the minimum value of the coefficients of the fourth 


order characteristic equation are defined as follows: 


2 
T 


3 


s" (2t. Ja, s +(1+a, o "s^ (2t, a, )u Asta o t = 0 (F20) 


Similar examination of higher order characteristic 
equations shows the relationship between a's, © and On) 
and their coefficients. Appendix B lists the various forms 
of these coefficients, including minimum stability values, 


for model equations up to seventh order. 
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As & result of the considerations of static and dynamic 
specifications we can now write model characteristic equations 
for a system which will exhibit a specified response. Con- 
siderations of chapters II and III have given us an actual 
characteristic equation which must be forced into this 


Specified configuration by compensation. 
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CHAPTER 5 


OF THF CHARACTERISTIC EQUATION 


5.1 Introduction. In Chapter 3 we formed the characteristic 
equation which represents the components in our servomechanism. 
The static and dynamic specifications were related to coeffi- 
cients of model characteristic equations, and these are tabu- 
lated in Appendix B. The third phase of our compensation 
procedure consists of determining the effect of compensation 
on the coefficients of the characteristic equation. This effect 
can be considered under two headings: 

a) The Effect of Compensator Location 

b) The Effect of Compensator Type. 

The purpose of this chapter is to show how to quickly deter- 
mine each of these effects. The relationship between them 
and the coeffitients will be examined. A notation will be 
introduced which permits rapid evaluation of the effect of 
all possible compensators, without repeated expansions of the 
characteristic determinant. 

The separation of the two effects is possible largely 
because of the unique manner in which the addition of a com- 
pensator alters the characteristic determinant. Before con- 
sidering the two effects separately we will examine the effect 


of compensation on the determinant of a general, four-node 
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uncompensated servomechanism. 


5.2 The Effect of Compensation on the Characteristic Deter- 
minant. In studying the effect of a compensator on the charac- 


. teristic determinant, we will use the four node system in 


fipure 5.1. 





Figure 5-1 
The characteristic determinant for this system is: 


L O 0 Gna 


aD ` Da 
-G A 1 O 0 
ZN = 
0 0 “Gac a 


The characteristic determinant is single valued, and so its 
value is independent of the method used in expanding it. 


For example: When expanded by cofactors of the fourth column 


EMO CEu Jm i o Eo] 
O -GB l -G * 
0 O “Gio O -G 





2d 


A = GapGpa Ga app Ge c*i 


A similar expansion by cofactors of the third row yields: 


l O Gap? pa l 0 Gap "pa 


a > E TO -GA O O d elas as I 0 


0 -G al 0 0 1 


eb 


A = GypSanSnalechan * 1 


Note that both expansions give the same result, and indeed, 
any expansion would give this result. 

If we were to add a compensator to any path previously 
open, we would put a transfer function in a determinant posi- 
tion previously zero. As an example, consider closing the 
path from the output of block Cop to node "a" with an element 
whose transfer function is Gha’ The characteristic deter- 


minant of this modified system is: 


l1 GaG 0 G.A e 


bB"Ba dD"Da 

A, = a ji 0 0 
0 -G op 1 0 

0 0 “Gac | 


Expanding by cofactors of the first row: 
1 O O -G 0 O 


eee) 1 0 so a E 





+0 


“Coa 1 0 
(GapGpa? O0 — -Gpp 1 
0 0 -G c 


Am = 1 * Sypfpafan + GapGpaSanSpp foc 
This is merely the old /\ with the term Gy gGp4 04, added. 


A similar effect is observed is an existing path is 
paralleled. Suppose an element whose performance function 
is GAA was put in parallel with Goat This would make the 
transfer function between nodes "a" and "b" equal (Goa t E 
Since no new position in the determinant was filled, the 


expansion remains: 
- ' 
An = (Gaa * Caa? Cbpecêan pa + 1 


A, = Saa% pp ecfapfpa * 1 * Gia pp ecO appa 


This again is the old /N with a new term added. We can gene- 


ralize the above observations as follows: 
A, = A+ G, Ay; (5.1) 


En TOW, qun column of 


where Gij is the element added to the i 
the determinant as a result of adding the compensator. Aij 
is the signed minor (cofactor) of this determinant position. 
It is the effect of G; As; on the coefficients of the charac- 


teristic equation that we must determine. 


All the terms in equation (5.1) contain ratios of 
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polynomials in "s". These polynomials come from the products 
of terms Gun? Gym? etc. in the characteristic determinant. 
These transfer functions represent components of the servo- 
mechanism. By their physical nature, components normally 
contribute a time lag to the system. This lag appears in 

the transfer function as an (sta) factor in the denominator. 
The expansions of /\ and a consist of unity plus a series 
of products of such transfer functions. Because of the pre- 
valence of lagging components these products will have nume- 
rators of order less than (or at best equal to) their denom- 
inators. When the terms comprising ZN are placed over a common 
denominator, the numerator so formed is the usual form of 

the characteristic equation. The common denominator will be 
the product of all the non-repeated (sta) factors, and the 
highest power of s in any of the individual product's denom- 
inators. The numerator contains unity times this common 
denominator. It is this contribution to the numerator that 
contains the highest power of s, and thus sets the order of 
the characteristic equation. 


With these thoughts in mind, consider again equation (5.1): 


Ay T A * GLAS (5.1) 





Modified Char.Eqn. . Char.Egn. 
i Zorn, + 6,44; (5.2) 


[C:P:m] Char. Fan. +G 
C.D.A 


Modified Char. Eqn. 


131; PAn 
(5.3) 


From this we see two types of changes to the coefficients of 
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. the characteristic equation: 


a) If pan is not unity, i.e., if the compensator 
added new (sta) factors or a higher power of s to the common 
denominator, the old characteristic equation will be multi- 
plied by the product of these added terms. This will raise 
the order of the characteristic equation. 

b) The coefficients of s* in Gy AAs 4D, Will be 
| added to the coefficients of s in [Dan | characteristic 
equation. 2 

Had the purpose of this chapter been merely to show the 
effect of a compensator on the coefficients of the charac- 
teristic equation, it would be fulfilled. One would only 
have to apply equation (5.3) for each compensator-path 
combination. For an "n" node system and "m" possible com- 
pensators, this requires (nx n x m) applications of equation 
(5.3). This is a process that FINES many, many hours of 
algebra. It usually just isn't done. The designer relies 
on his previous experience to help him pick a usable path and 
is likely to overlook one which would have done the job 
easier. The remainder of this chapter will develop a tech- 
nique for examining all possible paths in a reasonable time, 
and for displaying this information in a useful manner. 

5.3 The Effect of Compensator Location. We have seen 

(equation 5.1) that adding a compensator has the effect of 
adding a tern, Gj jÂ? to the characteristic determinant. 
Factoring the performance function of the compensator from 


the term Gs As; will leave the effect of closing the path 
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in which the compensator was inserted. It is this effect 
that we call the "Effect of Compensator Location". GrG" 
EG 
It is the term which would appear as characteristic deter- 
minant element i-j, if the path represented by that element 


were merely closed. 


Gay = Gyn ties, = 


Rewriting equation (5.3) for this condition: 


_|CDAm 
Modified Char. Eqn. -| Char. Eqn. + GEA ¿DA my 


as 
(5.1) 


Since there is no denominator in the performance function 
Go = 1, there are no new factors added to the common denomi- 


nator by the inclusion of GEA, y" Dan’ PA » 


Modified Char. Eqn. = Char. Eqn. + Ay ¿Da (5.5) 


Medified Char. Fan. - Chars Eqn. = Gr 40DA (5.6) 


(Modified Char. Eqn. - Char. Fqn.) will be a polynomial in s 
containing only those powers of s whose coefficients were 
changed by closing the path. If we could quickly determine 
the powers of s in GA; jA » for all possible paths, we 
would know what coefficients of the characteristic equation 
were altered by closing those paths. 

To accomplish this we will introduce a new notation. 
This notation will describe the powers of s contained in the 


polynomials which comprise the terms in equation (5.6). 
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Since these polynomials in s arise from the multiplication 
or division of factors (sta) or "s", they will contain all 
the successive powers of s between the highest and lowest 
power present. The symbol so used will be called an "Order 
Number", and will consist of two numbers, one Arabic and 
one Roman, i.e. (A,R). The highest power of s in the poly- 
nomial is and the lowest power of s is BS. All powers 
of s between these are present. For example: 
If (Order of P(s)) = (3,1), then P(s) = a, 5 ta s ras s^ ea, s. 

The algebra of these Order Numbers is much like that of 
exponents. When two polynomials are multiplied, their Order 
Numbers add to give the Order Number of the product. If 
one polynomial in s is divided by another, the difference 
in their Order Numbers is the Order Number of the quotient. 
Order Numbers add and subtract as follows: 

(A4 ,R4) t (45,R5) 7 (A,tA5, RYtR,) 
As an example of this use of Order Numbers, consider: 

(sta) (s“+bs) = 5° + (a*b)s^*abs 

Order Number of (sta) = (1,0) 

Order Number of (s^*bs) - (1,1) 

er + (IST) - (2,1) 

Order Number of s? t (a*b)s^*abs saa 1) 
Since the terms of equation (5.6) are polynomials in s of 
the type representable by Order Numbers we can rewrite this 


equation as follows: 


(Coef. Indicator) = (Order of G,) + (Order of Ay y? + 


(Order of DA? (5:7) 
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Before examining the use of equation (5.7), it is well 
to become more familiar with these Order Numbers, and to see 
how they are found for the terms of the equation. 

The Arabic number equals the highest power of s arising 
from the combination of the terms of the form (sta). 


Example 1: In (sta)(stb)(stc), the Arabic number is 3. 


Example 2: In CIEL the Arabic number is 1. 
stc 


The Roman number equals the net power of s arising frön 


terms of the form sa, 


3 


Example 3: In s^, the Roman number is III. 


Example 4: In mo the Roman number is (m-n). 
S 


As examples of polynomials in which both Arabic and Roman 
numbers occur: 


Example 5: The order of (s*a) is (0,I) 
s(s*b) 


Example 6: The order of s^ (s*a)(s*b)(s*c) ISO LPS 
Example 7: The order of O is (-m,-a@) To simplify the 
notation, this will be designated (-,-). Note that (A,R)+ 
(-,-) = (-,-). Any sum containing (-,-) has no value. It 
is as if you multiplied by zero. 

The Order Numbers comprising the left side of equation 
(5.7) are obtained as follows: 

(Order of Gr): In determinant positions corresponding 
to feedback paths, Gr, is the transfer function of the block 
immediately preceding the feedback pickoff. It would be of 
the form of Example 5 above. In determinant positions corres- 


ponding to feed forward, or parallel to a direct path, it 





46 


equals unity. The Order Number of a constant is (0,0). 

(Order of As 3): This could be found by expanding the 
minor, Ars and then applying the rules exemplified above. 
A simpler method consists of replacing each existing term 
in the characteristic determinant with its Order Number. 

The minors are then expanded by the usual determinant rules 
with the following exceptions: 

1) The Order Numbers are added where the terms they 
represent would have been multiplied. 

2) Arabic and Roman numbers are added separately. 

3) Ihe largest Arabic and the largest Roman numbers 
resulting from these additions comprise the Order Number of 
the minor. 

(Order of CDA): The common denominator of /\ will have 
already been found in the determination of the characteristic 
equation (Chap. 3). Apply the rules above to this denomi- 
nator as in Example 6. 

With the use of the Order Number technique, a matrix, 
showing the results of equation (5.7) for each path, can 
be formed. This “Location Matrix" is constructed as follows: 

a) Replace the existing transfer functions in the Char- 
acteristic determinant by their Order numbers. 

b) Evaluate the (Order of A; y? for each i-j position 
usable for compensation. (Normally all positions) 

c) To this (Order of Ay y) add (Order of G,) for the 
position concerned. 

d) Add (Order of Gi/A,,) to (Order of CDA) and display 


the resulting Order Numbers in a matrix. Each Order Number 


E tu — 
^. >> $ 
EJ 
n BEEN mE 
Y ==> 


P ES y 
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will occupy a position in the matrix corresponding to the 
position of its path in the Characteristic Determinant. 
This matrix will indicate which coefficients can be 
changed by closing a given path. 
To clarify the ideas presented in this section, the 
location matrix for a sample problem will be formed.  Con- 


sider the system of Figure 5-2. 





Figure 5-2 


This system has the following characteristic determinant: 


1 0 0 E, 
S ( sta ) 
-K 1 0 0 
2 
-(stc 
P ES 1 s 
0 0 ES 7 
(s+b) 


Expanding this in terms of elements of the fourth column: 
Kj KK, ( stc ) 


A= 1+* seta) (stb) (sa) 


A = s(sta)(stb)(std)*K K-K, (stc) 


s(s+ta)(s+b)(s+a) 
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Expanding the numerator and introducing standard symbols for 


coefficients: 
n s d “+a Sta 
Char.Ean. - 2 I 0 
MN s(sta)(stb)(std 


Applying the rules for formation of a Location Matrix: 


Step a: (0,0) (-,-) (-,-) (-1,-1) 
A ha (0,0) (0,0) (-,-) (-,-) 
(-,-) (0,0) (0,0) (7, -) 


(-,-) (-,-) (150) (0,0) 


Step b: A}; ~- (0,0) Ayo - (0,0) Arz - (0500.78. - E1507 
As E (-2,-1) Apo F (0,0) AVE = (0,0) As PX (-1,0) 
Ary 7 25-1) As, - (-2,-D) Azz - (0,0) Aa, - (71,0) 


Ayı P (-1,-I) Ayo = (-1,-1) Ay 3 ES (-1,-1) Zins = (05:09 


The Order numbers of the Gr's are as follows: 


G G G G 

258 (00,0) ^1» - deko) 13. > (1,0) u = Ce 
Go, G G G | 

Pai - (0,00 le2 - (0,0) %23 - (-1,0) “ak - (-1,-1) 
G G G G 

A (sc) 532 > (000) "33 - (21,0) m Ca) 
G 


G G G 
= (o) ge (0,0) Au - Ce 


Step e: (Order of GA; y? - position in the array denotes 


subscripts. 


(0,0) (0,0) (71,0) (-2,-1) 
(-2,-1) (0,0) (-1,0) (-2,-1) 
(-2,-1) (-2,-1) (-1,0) (-2,-1) 
(-1,-1) (-1,-1) (15-1) (-1,-I) 


(Order of CD A) equals (3,1). 
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Bien [3,1 ED a) (1,0) 
Gu) (3,1) (2,1) (1,0) 
Location Matrix - (1,0) (1,0) (2,1) (1,0) 
(220) (2,0) (2,0) (250) 


Step b, finding the Order of each minor, merely entails 
scanning the auxiliary determinant of step a to find the 
largest Arabic and the largest Roman sums present. The 
presence of numerous (-,-) terms in an uncompensated system 
makes this a few minutes work. This is in comparison to 


the hours of algebra required to evaluate G4 444 40D, for 


every possible compensator. 

The Location Matrix is the sum of the terms of Equation 
(5.7) for all possible compensation paths. The information 
in the Location Matrix, modified by the effect of compensator 
type, will indicate changeable coefficients of the Charac- 
teristic Equation. 

93. The Effect of Compensator Type 
5.3.1 Introduction 

Examination of the type of compensator yields informa- 
tion on: 

a. Which coefficients of the Characteristic Equation 

will be changed, and 

b. the extent of that change. 

It is the Order Number of the compensator in conjunction 
with the Location Matrix that indicates which coefficients 
will be changed. The performance function indicates the 


magnitude of the change. 
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Table V-1 lists the compensator types to be considered, 


with their Order Number and Performance Function. 


Table V-1 Elementary Compensators 


PERFORMANCE 
ORDER NUMBER FUNCTION 


K 









TYPE 






Sensitivity 





S 






Differentiator 





First Order Lead Filter (sta) 













Second Order Lead Filter (s**bs*c) 
Integrator 1/s 
First Order Lag Filter 1/(sta) 






1/(s*+bste) 





Second Order Lag Filter 


5.3.2 Determination of Which Coefficients are Changed 


The Order Number of the compensator, when added to the 
element of the Location Matrix corresponding to its path, 
gives a new Coefficient Indicator. As previously defined, 
the (A,R) pair resulting from this addition indicates which 


coefficients of the Characteristic Equation are modified; 


i.e., SACR? 3A-14R? =; ag’ 
There are three situations in which addition of a com- 
pensator will raise the order of the Characteristic Equation. 

If the order is raised, there will be two changes to the 
coefficients. One change will be to all coefficients due 
to merely raising the order. The change dictated by the 
Coefficient Indicator will be an additional change to the 

coefficients of the higher order equation. 


Before examining the quantitative effects of compensation 
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on the coefficients, the three cases where the order of the 

Characteristic Equation is raised will be considered: 
Case I - This situation arises when integrator com- 
pensation is used. If the addition of a con- 


pensator results in one of the terms of 


Loe too b> + Sap ie + - = = 
r T ro r 3 B 


having more integrations than the overall 
uncompensated closed loop performance func- 
tion, the order of the Characteristic Equation 
will be raised. In the above expression, In 
signifies loop performance function. TT 


signifies products of loop performance 
functions taken k at a time. That the order 
of the Characteristic Equation is raised in 
this situation is apparent if the Character- 
istic Equation is derived from the system per- 
formance function obtained by application of 
Mason's Roe 


A simple example of this effect may be 


seen by considering the following system: 


Block A Block B 
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Char. Egn. . S * "re 
EM S 
K,K, = O 


Char. Hén. = 5 + “12 
For this system the Location Matrix is: 


(0,1) (0,0) 
L. M. - 
(0,0) (0,0) 


Consider integration feedback from the output of 
Block B to Node b: 


Block A Block B 








-K l+ 2 

l 2 

S 

A= 1 2.03 








23 


2 


Char. Fan. _ s + er > is K2 
za S 


— = 
Char. Eqn. = s + ESE + K, -0 


Notice that the order of the Characteristic 
Equation has been raised by this combination 
of compensator and path. If the integration 
had been fed back from the output of Block A 
to Node a, a loop performance function higher 
than first order would not have been generated 
and the order of the equation would not have 


been raised. In this case, 





Char. Fan. _ s + (Kj) (1+K>) 
As — S 


s + K, (1*K5) = © 


Char. Eqn. 


Table V-1 gives an Order Number of (0,-1) 
for this compensator. Adding this to the 
proper element of the Location Matrix gives: 


Fed from B to b: (Coef. Ind.) = (0,-I) 


Fed from A to a: (Coef. Ind.) = (0,0) 
Since the Characteristic Equation contains no 
negative powers of s, the first indicator (for 
B to b) appears to be contradictory. It is 
this very contradiction that yields the 


information on equation order. The order of 


the Characteristic Equation must be raised 





Rudié-I - 


: Case II - 
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by the amount that the "R" of the Coefficient 
Indicator is negative. That this occurred is 
borne out by comparison of the Characteristic 
Equations in the above example before and after 
this compensation. In the case of feedback from 
A to a, the "R" of the Coefficient Indicator 
accurately reflects that there is no change in 
order. 

The order of the Characteristic Equation is 
raised by the amount that the Roman Numeral of the 
Coefficient Indicator is negative. 

Another order raising situation arises when 


any of the terms of 


IE E A EL t ERL +- 
p E r2 T DES D 


when placed over a lowest common denominator, has 

a numerator of higher order than the uncompensated 
equation. Again, this is apparent if the Charac- 

teristic Equation is derived from the system per- 

formance function obtained by application of 


Mason's Rule. 


For example, if the sample system is compen- 


sated as follows: 


Block A Block B 








p5 


thems /\ =| 1 Ko 
S 


-K, —1* Ko(s^ bsc) 
S E 


AA + BR + Ko (s^*bs*c) 
S S 


Char. Ean. _ Kos“+ (obs t (Kee + KIKO 
AD S 


2 
Ghar. Eqn; = s^ + (Kabri } + (K,+c) = 0 
2 


Adding the Order Number of this compensator, (2,0), 
to the proper Location Matrix element gives a 
Coefficient Indicator of (2,0). Again there is an 
apparent contradiction in that the highest coeffi- 
cient indicated to be changed is not present in 
the uncompensated Characteristic Fquation. This 
contradiction indicates that the order of the Char- 
acteristic Equation must be raised. 

Putting this compensator in the path from the 
output of Block A to Node a yields: 


N = Char. Fan. _ Ky s? + Ky bs“ ENEE Ls IO 
SN S 
Char. Eqn. = K ED dE bs* + (K,c+1)s + K,K, = O 
e. qn. i di 2l 1^2 


The Coefficient Indicator in thís case is (2,1), 
which again indicates a raise in the order of the 
Characteristic Equation, the coefficients to be 


changed being 845 85 and az» 





56 


Rule II - The order of the Characteristic Equation is 


raised by the amount 
(Auer. mua. AD of CDA 


Note: If an occasion arises where both Rules I 
and II apply, first raise the Order of A by 
the amount indicated in Rule I, put "R" of the 
Coefficient Indicator to zero, then, if necessary, 
apply Rule II. 

Case III - A lag filter will always introduce its denom- 
inator in EUN . This will raise the order of the 
common denominator by an amount equal to the order 
of the filter. Since in practical servos, the 
order of the Characteristic Equation equals the 
order of the common denominator of A, it too 
will be raised by this amount. 

Rule III - The order of the Characteristic Fquation is 
raised by the amount the Arabic Numeral of (Order 


of the Compensator) is negative. 


5.3.3 Determination of Quantitative Changes to Coefficients 
Me Modified Characteristic Equation, Equation (5.3), 


consists of two terms: 


Mod. Char. Fan. = Am (Char. Fan.) + ¡DA 
DA, J 
This indicates that there are two contributions to the change 
in the coefficients. The first term's effect is predictable 


in the case of simple compensators. It is the amount all 
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coefficients are changed when the order of the Characteristic 


Fquation is raised by applying Rules I and III above. 


In Case I, DAn t sIRl where the "R" 
OL IN d is the Roman 
Numeral of the 
Coefficient 


Indicator. 


In Case III, ÜD/Am | (s*a) or (s^*bs*c) 
RAN 

depending on the 
order of the fil- 


ter used. 


In these cases, the initial effect on coefficients can be 
readily calculated. 

Ihe change to coefficients BAR? SA-14R? etc., as indi- 
cated by the Coefficient Indicator can only be found by eval- 
uating GA DA for the path concerned. Evaluating the 


change on coefficients caused by raising the order of the 
Characteristic Equation in accordance with Rule II can only 


be determined by evaluating G5 444 PAn" 
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CHAPTFR 6 


SELECTION OF COMPENSATOR 
l1. Introduction 


The problem to be considered in this chapter is that of 
choosing the combination of compensator and location that can 
most easily change the characteristic equation of the system 
to one established by the specifications. This requires a 
comparison of the coefficients of the uncompensated charac- 
teristic equation to those of a minimum model equation., The 
actual coefficients which exceed a minimum allowable value 
will be used to select a working model equation. The work- 
ing model equation permits calculation of a compensated char- 
acteristic equation consistent with the specifications. 

There are four equations of interest in this chapter. 
Ihe uncompensated characteristic equation describes the 
servomechanism before compensation. Its coefficients are 
set by the basic hardware. It is these coefficients which 
must be changed. The minimum model equation is the charac- 
teristic equation of a system with a transient response which 
has: 

a. A specified dominant damped oscillatory mode, 

b. A minimum number of first order modes, (none or one), 

GN Non-dominant undamped oscillatory modes. 


No first order modes will exist for even ordered systems. No 
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non-dominant oscillatory mode will exist for a third order 
system. Its coefficients are those of the characteristic 
equation of a system verging on instability. The working 
model equation is the characteristic equation of a system 
with a transient response which has a specified dominant 
damped oscillatory mode. Its non-dominant modes will be 
damped oscillatory and/or first order. Its coefficients are 
those of a characteristic equation of a stable system. The 
compensated characteristic equation describes the servo- 
mechanism after compensation. Depending on its stage of 
development, its coefficients may be purely numerical or 
may contain unspecified parameters. Compensation which 
changes the coefficients of the uncompensated equations to 
those of the compensated equation will result in a system 


with the specified response. 


2. Choosing the Minimum Model Fquation 
There are three factors which determine the choice of 

a minimum model equation. They are: 
dug It must be of the same order as the uncompensated 
characteristic equation. 
er The quadratic made up of the specified Z and w 
must be a factor of the model. 
3. The dominance factors, a's, must be sufficiently 
large to insure that the response of the model has the 
desired correspondence with the response of the quadratic. 
Ihe minimum model equations in Appendix B are so con- 
structed that substitution of the desired č, œw. and a's 


n 
will give them the properties listed above. 
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3. Determination of the Compensated Characteristic Equation 


The compensated characteristic equation must have coef- 
ficients at least as large as those of the minimum model. 

If its coefficients are less than those of the minimum model, 
the system must be unstable. A comparison of the uncompen- 
sated characteristic equation to the minimum model equation 
will show which coefficients must be changed. There are 

two possible conditions arising from this comparison. Either 
all of the coefficients are below the minimum, or a number 

of them are already sufficiently large. In the first case 
substitution of non-dominant damping ratios and/or acceptable 
dominance factors into the working models of the proper order, 
as given in Appendix B, will set the compensated characteristic 
equation. The ratio of a, to ag in the compensated charac- 
teristic equation may be set by static specifications. This 
may limit the choice of a's and non-dominant C's. 

In the case where some or all of the coefficients are 
already above the minimum, several factors must be considered. 
These ares 

l. The interdependence of the coefficients of the model 
equation. 

2. The position of the satisfactory coefficients. 

3. The effect of static specifications. 

The interdependence of coefficients is expressed in the 
working model equations of Appendix B. In an nos order model 
the n coefficients are functions of the specified ¢ and Qn 
and of (n-2) dominance factors and/or non-dominant C's. 


Ihis gives rise to n equations in n unknowns. Two of these 
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unknowns, & and a, of the dominant quadratic, are set by 
specifications. The remaining (n-2) degrees of freedom 

may be constrained by fixing either coefficients or a's. 

If any (n-2) coefficients of the uncompensated characteristic 
equation are above the minimum values, they will determine 
the required value of the remaining two coefficients. If 
fewer than (n-2) coefficients are satisfactory, the 
remaining degrees of freedom may be fixed by an arbitrary 
choice of a's. Required values for the remaining coef- 
ficients can then be calculated. These satisfactory and 
calculated coefficients comprise the compensated character- 
istic equation which will require minimum change in coeffi- 
cients. The satisfactory coefficients should first be com- 
pared to the model which contains all real non-dominant roots. 
If they are not sufficiently large to fit this model, then 
the oscillatory working model equation must be used. 

The second consideration arises from the fact that all 
compensators affect successive coefficients in the character- 
istic equation. As developed in Chapter 5, the independent 
coefficient changes are found by evaluating G3 444 4C DAR 


This will be in the form of Ks'"(s*a)(s*b) etc. When 
expanded, it will contain successive powers of s. This fact 
is also revealed by the Coefficient Indicator. Because of 
this fact, the satisfactory coefficients used to determine 
the compensated characteristic equation should be chosen so 
that the coefficients remaining to be changed are successive. 
The static specifications placed lower limits on some 


coefficients of the compensated characteristic equation. 
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These limits may force some coefficients to be larger than 

the dominance factors dictate. The coefficients so affected 
must be included in the (n-2) degree-of-freedom considerations 
above. 

Comparison of the desired Coefficient Indicator and the 
Location Matrix will determine the type of compensator required 
to change the necessary coefficients. If it raises the order 
of the equation, the abpve procedure must be repeated start- 
ing with the minimum model equation of the higher order. The 
uncompensated equation to be compared with the model can be 
formed knowing the order raising effects described in 


Chapter 54 


4, Evaluating the Parameters of the Compensator 


The Gs As; of the required compensator is added to the 
old characteristic determinant. This results in a charac- 
teristic equation in which the coefficients requiring modifi- 
cation contain the parameters of the compensator.  Equating 
these coefficients to those of the compensated characteristic 
equation will determine the value of the parameters. The 
compensator must have sufficient degrees of freedom to 


match the independent coefficient changes required. 
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CHAPTER 7 


COMPENSATION DESIGN PROCEDURE 


7.1 Introduction 


Ihe advantages of the technique described in the previous 
chapters are three. First, it involves no graphic presenta- 
tions as required by other methods. Second, it does not depend 
on the intuition of the designer for its choice of compensation 
to achieve the desired result, Finally, it requires little 
of the extensive algebra inherent in other analytical methods. 
Since there is no intuition or experimental graphing involved, 
an ordered analytical approach to any compensation problem 
is possible. 

This method divides naturally into four phases, which 
are the subjects of Chapters 3 through 6. These chapter 
divisions were made primarily in order to explain the tech- 
nique. They also follow in the proper logical order for 
application to the solution of a problem. The step-by-step 
procedure in this chapter thus serves as a summary of the 


paper as well as a guide for the compensation designer. 


7.2 Step-by-step Guide to Compensation 


Given: I. A servomechanism for which the direct path 
components have been selected. This servo may have 
been partially compensated or may just have unity 


feedback from output to input. 
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II. Static and Dynamic Performance Specifications 
to be met by the system. 
Problem: To design compensation to meet the specifications. 
Solution: 
Phase I: Forming the characteristic equation. 
A. Arrange the system in a standard block diagram. 
COMMENT: Combine the components between usable 
pickoffs and nodes into single blocks. 
B. Form the characteristic determinant. 
C. Write the characteristic equation by setting 
the expanded characteristic determinant equal to 
Zero. 
COMMENT: Note the Order Number of DR 
Phase II: Relating Specifications to Coefficients of 
Model Characteristic Equations. 
A. Express static specifications in terms of 
restrictions on ao and ay of the compensated 
characteristic equation. 
B. Interpret dynamic specifications in terms of 
second order response characteristics; i.e., © and 


ei? 


C. Choose minimum acceptable dominance factors 
(a's and non-dominant C's) to insure second order 
dominance. 
Phase III: Determining the Effect of Possible Compen- 
sators on Coefficients of the Uncompensated Character- 


istic Fquation; i.e., Forming the Location Matrix. 
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COMMENT: Chapter 5 shows the details of this phase. 
Phase IV: Choosing the Compensator. 

A. Form the minimum model equation using a's, €'s and 

c, found in Phase II. | 

COMMENT: All coefficients of the compensated equation 

must be greater than those of the minimum model to 

insure stability. 

B. Choose the Working Model Equation. 

l. Insure compliance with dynamic specifications. 
Put sufficiently large a's (Phase II-C) and the required 
C and aw, (Phase II-B) in the working model equations of 
the proper order. 

2. Insure compliance with static specifications. 
Compare 81/86 ratio, of the various forms of the working 
model equation. Choose for a working model that form 
which, with allowable non-dominant zetas, (Phase II-C) 
has the smallest coefficients that meet the specifications. 

3. Remove the effect of static specifications from 
working model equation. In the event that the minimum 
model equation will not meet static specifications, the 
basic system configuration must be changed by addition 
of a component which reduces following error to zero. 
This will invariably mean the addition of an integration. 
Compensation of this type will introduce a zero in the 
system performance function. Set this zero to a value 
which will remove its effect on the dominant mode of the 
response (See Figure A-2). The design procedure must 


be started again with the revised system. 
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Form the compensated Characteristic equation. 

1. Compare the uncompensated characteristic equa- 
tion to the working model. 

2. If n-m-2 coefficients of the uncompensated equa- 
tion are above the working model coefficients and 
leave successive coefficients to be altered, use 
them to fill some of the n-2 degrees of freedom. 
Set m arbitrary a's or non-dominant C's. Compute 
the m*t2 remaining coefficients of the compensated 
characteristic equation. If no coefficients are 
above those of the working model, it is the com- 
pensated characteristic equation. 

Select the Compensator. 

1. Form the desired Coefficient Indicator from the 
coefficients in C-2 above which must be changed. 

2. Compare this Coefficient Indicator to the Loca- 
tion Matrix to find the simplest possible compen- 
sator. The compensator selected must introduce 
enough additional variable parameters so that the 
total number equals the number of coefficients to 
be changed plus the number of zeros of the system 


to be set (See COMMENT) 


3. Form EE , which is the system perfor- 
In 


mance function relating the actuating input at 


node "n" to the putput of block G This 


jJ" 
expression will contain the adjustable parameters 


UW. 


of the compensator. The compensated characteristic 
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equation will be the denominator of this function. 
Equate the parametric form of the coefficients of 
the denominator of the system function to the coef- 
ficients of the compensated characteristic equation. 
Also set the parameters of any numerator zeros to 

a value which gives a satisfactory zero (i.e., 


sufficiently large a). 


COMMENT: If the only possible compensator raises the 


order of the characteristic equation, apply the 
order raising rules of Chapter 5 to the uncompen- 
sated characteristic equation. This order raising 
will put a zero in the system function. Repeat 


phase III using higher order model. 
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CHAPTER 8 
SAMPLE PROBLEM 


8.1 Introduction 

The problem to be worked in this chapter is a modification 
of an analysis problem presented by Chestnut and Mayer. 1? Ihe 
problem, as presented in Chapter 13 of the reference, consisted 
of determining the performance of a position control system 
subjected to a load disturbance. Compensation had already 


been accomplished. It consisted of feedback of the first and 


second derivatives of the output quantity. The signal fed 


back was modified by a filter of the form Kae sta e 
(stb) (ste) (std) 


The specifications were estimated from the performance of the 
system as compensated in the reference. 

8.2 Statement of the Problem 

Given: I The system to be compensated has the primary require- 
ment of position control. However, it must respond properly; 
within specified limits, when subjected to a designated 
disturbance load torque. The main path components consist 


of the following: 





Component 
Preamplifier 


Amplifier, Amplidyne 


Motor and Load 


Disturbance Load/Voltage 
Transducer 
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Transfer Function 


Ki not specified 


PIE 


ka = 6855 volt/ft.1b. 


T 
II Static Specifications 
Steady State Error in response to step = 0.0 
Following Error in response to ramp = 10^ rad/rad/sec 


Maximum output in response to load torque - 10 rad/ ft, 1p. 


Dynamic Specifications: 


Peak Overshoot 


Settling Time to A of steady state 


in response to step 


= .6 sec 


Problem: Design compensation such that the system will 


conform to specifications. 


8.3 Solution of the Problem 


Phase I: Forming the characteristic equation. 


A. 





Load Torque 
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B. 
1 0 0 4.78 
s(s*9.1) 
NE -K, 1 0 0 
0 -1500 p 0 
2+31.5 
O 0 -050 dE 
s+t17. 
Ce 6 


a 6.09 x 10 K 
os Lot se E T 
6 


A=8 Stil. £t1716)(8$*9,1)t- (6409 "* 10") K 
s(s+31.5)(s+17.6)(s+9.1 


A= s^ 58.2531 0017 «506084 (6.09 EXE K] 
s(st31.5)(s*17.6)(s*9.1 
y à 2 6 * 
Char. Hon: © + 58.267 + 10008 + 5040s + (6.09 x 10 IK] z 0 
(Order of CD A) = (3,1) 
Phase II: Relating Specifications to Coefficients of Model 
Characteristic Equation. 


A. Static Specifications 


l. Steady State Position Error in response to step = O, 


a, = Limit [Sa A DA] 


S0 
o = Limit[_4.78 (1.275 x 107) : 
S 5 O[s(s*9.1) (s*31.5)(s*17. 


(549.1) (431.5) (5917.6) 


w 
i 


a 


6 
5 (6.09 x 10 IK] 


The value of ao is independent of this specification 
because of the unity feedback. 


2. Following Error in response to ramp = 107+ sec, 
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Error. = 107 ——2040 __ 


(6.09 x 10 IK] 


K, > —dÀ9:0 = 8.28 or = < 107* 
6.09 x 10 O 


3. Response to Load Torque = 107+ rad/ft.lb. 


k > Daun! Giray CDA 
ao 


1075 » Limit | 855) (4.78) (s)(s+17.6) (s+31.5) (s+ aj 
as, e © (6.09 x 10?)(s)(s*9.1)K, 


107" > (4827) 06.78) 17.62 (31.2) 
(6.09 x 10 IK] 
Ki > 3.43 
B. Dynamic Specifications. 
From reference 13, Peak Overshoot of 1.2 corresponds 
to (.6 > % > .5). We will choose the conservative 
č = 66. 
From reference 10, Error is reduced to within ICH 


o 
of final steady state value in four time constants, and 


Settling Time = + = .6 sec. 
l CO 


w = = = 11.1 rad/sec 


Co Choosing Dominance Factors. 

From Appendix A, 4, 2 3 and a, 2 2 are acceptable 
values of dominance factors. If non-dominant roots are 
complex, use a > 3 and t5 2 3. 


Phase III: Determining the Effect of Possible Compensators. 
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CORO) (-,-) (-,-) (-1,-I) 
(-,-) (0,0) (-,-) (-,-) 
Auxiliary A = 
(-,-) (2190) (050) (-,-) 
(7,7) (-,-) (-1,0) (0,0) 
(0,0) (0,0) (-1,0) (-2,0) 
(-39=1) (0,0) (-1,0) (-2,0) 
Order Numbers 
= „ee (0,0) ( 29985 
of Aij's 
(-1,-1) (-1,-1) (-2,-I) (0,0) 
(0,0) (-1,0) (-1,0) (21,1) 
(0,0) (21,0) (71,0) (-1,-1) 
Order nee. (0,0) (0,0) (-1,0) (-1,-I) 
"t hs (0,0) (0,0) (0,0) (1, 


(Order of CD A) = (3,1) 


BI (AND (dab) (0,0) 

(0,0) (257) (I) (0,0) 
Location Matrix: 

(ds) (180) (2,1) (140) 

(240) (2,0) (1,0) (2,0) 


Phase IV: Choosing the Compensator. 


E _ For a4- 
e CAL =e 


2 2 
a5) = (l+a, )&, = 13233 


A. 


2 u 
a, = 2%)%, @ = 14,080 
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Inn, 
ig.“ +, @ = 137,000 


aa + 133386? + 12234" * 1080s + 137000 = O 
Note that a, and 81 of the uncompensated system are 


below minimum values. 


B. aj for minimum model - .103 
ao 
aı for all complex non-dominant roots = .126 
ao 
aj for all real non-dominant roots - .156 
a 
O 


Since a/a, from the model equations can never be 
less than .103, the condition on a,/a, imposed by static 
specifications must be relieved by making the following 
error independent of this ratio. This effect can be 
achieved by compensating so as to have zero following 
error. This implies making the servo a velocity control 
system. However, it may not lose its position control 
ability. Paralleling the main path with another forward 
path containing an extra integration will result in both 
position and velocity control. This effect can also be 
obtained by placing the extra integration in a feed for- 
ward path from the first to second nodes thereby utilizing 
succeeding integrations in both forward paths simultaneous- 


ly. Using this compensation, the transfer function, Gan? 
becomes Ki(s * Ko/Kj) 
S 
In addition to raising the order of the characteristic 


equation, a zero of the system performance function will 


be introduced at s - -KJ/K, . Set this zero such that its 
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a > 2, Say K,/K] = 25. It is now necessary to recommence 


the solution. 


Phase 1: Forming the characteristic equation. 


is m ue 





B. 
4.78 

ji 0 0 SGH 
K 

2s =) 1 0 0 

O 
0 aS n 0 
_ 850 

5 e s+17.6 š 


" A s? + 58.2s* + 1001s> + 5040s” « 6.09 x 10 Kis + 6.09 x 10°K 


8 (8995 1)65«37,6) (&*31 75) 
Char.Eqn. - «58.25 ‘+1001 s°+5040s+6 .09 x 10°K, s+6 .09 x 10°K,=0 
(Order of CDA) = (3,11) 


Phase II: Relating Specifications to Coefficients of Model 
Characteristic Equations. 
A. Statis Specifications: 
l. Steady State Position Error in response to a 


step = 0 
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6 
a, = Limit [aAa] = 6.09 x 10 K, 


s 70 


As before, the value of a, is independent 


O 
of this specification. 


2o Following Error in response to a ramp - O 


Hérer - sat 2 [1 - 6 ING 
oe | =e 


3 Léset 5. - 58.2.) » 10075 + poles 
8 Char. Eqn. 


Therefore, the values of the coefficients are 
independent of this specification. 
3. Response to Load Torque < 107+ rad/ft.lb. 


d 
ior > Limit p 


Char.Eqn.. 


107* > Limit (.855)(5.78)(s)(s*17.6)(s*31.5) 


S > 0 Char. Eqn. 


107" > a 


Therefore, the specification on load torque 
disturbance is meto. 
B. Dynamic Specifications: 
Same č and wm, are required. 
C. The same dominance factors are acceptable. 
Phase IIIs Determining the Effect of Possible Compensators. 
The (Order of G,) for the G,, feedback position is changed 
to (1,-I) by the addition of the parallel integration. The 
(Order of CD A) is changed to (3,II). | 


(251) (2511) (1, II) (OT) 
New Location | 9:1) eM eM ee 
Matrix (ar) Gn) (2211) (15588) 


(251) (2,1) (1,1) (2,1) 
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Phase IV: Choosing the Compensator. 
A. Using i. 3 and a, = 2, the minimum fifth order 
model is: 
a? 2 senses” + 1527s? i l2, 100s* + 464,000s + 3,030, 000=0 
Since this is now a Type II system, it has zero 
error in response to inputs of position and velocity. 
The a,/a, ratio represents velocity error in response 
to an acceleration input. 
The requirement that the working model meet static 
Specifications is removed. 
B. 1. Since no working model had to be formed to con- 
sider the static specifications, the initial comparison 
will be made to the minimum model. With any easily 
realizable values for K, and K5, the ay, ay and a, of 
the uncompensated equation are above those of the minimum 
model. Since a, is the only fixed coefficient above 
the minimum model, it will be used as a satisfactory 
coefficient. K, and K, will be left unspecified. 
2. n-m-2-1 and m-5-2-1-2 arbitrary dominance 
factors. Solving the parametric form of the a, 
coefficient of the oscillatory working model, 
with 0473 and 0572, yields: 


(2%, * 2a * an) DSL 


1*2 
dan t 6%, + ee 
Lo = ‚534 

This non-dominant damping ratio is larger than needed, 


so Eo = 2.3 can be used to resolve for a larger PE 





n 


58.2 


a, = 2.24 


(1.2 +08 + a7)11.1 


Using these values for «a's and C's, the remain- 


ing coefficients can be found. 


us s 2820 
a, = 54600 
a, = 564000 
ig = 3415000 


The compensated characteristic equation is: 


E * 2820s? ar 546008" + 564000s + 3415000 = O 


de +58, 28 
C. 1. The coefficients to be changed are 845 85, 84 
and ao’ 84 and ao could be changed independently of 
additional compensation by adjusting Ky and Ko, but 
K,/K] has been set at 25. Therefore it is necessary 


that compensation provide for varying 84 and a, indepen- 


O 
dently of the other. In order to keep the coefficients 
to be changed suceessive, K, will be determined. From 
the characteristic equation developed in Phase I-C and 
the pica equation developed in Phase IV-C-2: 


6 


6.09 x 10° Ky = 3,415,000 
Therefore, Kk, 205725 = 10224 


The "semi-compensated" characteristic equation is now: 


eg E 58.25" + 1001s? ce 5040s" + 137,000s + 3,415,000=0 
The compensated characteristic equation is: 


y 


3? + 58.28 + 28205? + 546008” + 564,000s + 3,415,000=0 
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Comparison of these two equations gives the desired 
Coefficient Indicator to be (2,I). 

2. The compensator should be the simplest one 
which has sufficient adjustable parameters to change 
the required number of coefficients. The compensator 
must have three adjustable parameters since there are 
three coefficients to be raised. 

The compensator must be in a feedback path since 
feed forward compensators reduce the size of the coef- 
ficients they effect. Comparison of the desired Coef- 
ficient Indicator and the Location Matrix indicates 
that feedback position Gna requires the simplest com- 
pensator with necessary parameters, one with order 
number (2,0). This can be achieved with a transfer 
function of the form K(sta)(s+b) which will arise 


from the following: 





PF = KG I S SA.) 


3. a = A + CRAVE 


An 7 s?+58.25 '+100183+50408°+137,0005+3,415,000 


s (s+t9.1)(s+17.6) (s+31.5) 


a Q««78) (1500) (850) (are Aag Ac 
s(st921)(s*t17. 731. 
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G 65,4, = ME Y ease?) = -1290 ) 


Gra ap y | 137,000(s+2 
m s°+58.2s +(6.09x10 K,+1001)s° 





+(6.09x10°K, +5040) s°+(6 .09x10°K,) s+3 415,000 


Equating coefficients of the two forms of the compen- 


sated equations give the following relationships: 


a, = 6.09 x 10° K, + 1001 = 2820 or K3 = ¿000299 
a, = 6.09 x 10° K, + 5040 = 54,600 or K, = .00813 
a, = 6.09 x 10° Ke = 564,000 or Ke = .0925 


8.4 Checking the Solution 





The compensated system has the following block diagram: 





elt E. 
s(st9.1) 


LZ 
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ho 8 
l 0 0 s(s*9.1) 


" 4.78 -3 2 
long = UI p O | ceso S 


4,813x10 ^s*.0925) 


. .1590 
0 s T 0 
. 3D 
9 0 s+17.6 l 
Ly 


+28208>+54 ,6005°+564,0008 + 3.415.000 


A. = $7+58.2s 
" d (499.1) ($*3Y P). Eo T756) 


For an input of a unit step to node a and output from 


block Gap» the system performance function is: 


137,000(s+25) 


PF = „(8+214,9) (s°+205+1109) (s°+13.35+123.2) 


The response of the system in time domain can be found 
by the method of evaluating residues to be: 


-24 ,9t 


Ə ,(t) = 1+.0011e *.137e 719 * cos (31.8t-309.6) 


out 


41,h8e 70:96*5605(8.88t-1399.5) 


It is apparent that the contribution of the real root is 
negligible, being only on the order of Ps at its greatest. 
The contribution of the non-dominant oscillatory roots, though 
somewhat larger, is still less than one-tenth of the contri- 
bution of the dominant roots. It decays quickly and by t = .2 
adds only about ls to the value of the output. At t = .6 sec, 
the value of the output is about .98, a a deviation from 
final value. A1°/, deviation was specified. At t = .27, 
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peak overshoot occurs and equals 1.24; compared to the speci- 
fied 1.20. Failure to meet specifications exactly in these 
two instances is attributable to placing the zero too close 
to the origin. This indicates a need for further study of 
the effect of zeros upon the transient response, or for a 
more conservative choice of a's in connection with zeros. 

The solution obviously meets static specifications 
because of the introduction of the feed forward integration 
compensation. 

8.5 Summary of the Application of the Technique to this Problem 

The solution of the problem proceeded along straight- 
forward lines until it was discovered that a Position Control 
System could not provide acceptable following error. It was 
then necessary to stipulate compensation which would give the 
system Velocity Control characteristics. With this compen- 
sation inserted, the solution was then recommenced. The 
necessity of another compensator became apparent and its 
characteristics were determined. 

It is interesting to note the similarity between the 
compensation determined by this solution and that provided 
by Chestnut and layer. 1° Both involve feedback of first 
and second derivatives of the output. However, this solution 
replaces their elaborate filter, Ks“ (sta) y with 


(stb) (ste) (std) 


feed forward through the function and feedback of a signal 


< aly 


proportional to the output quantit 


o 
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CHAPTER 9 


CONCLUSIONS AND RECOMMENDATIONS 


9.1 Conclusions. Application of step-by-step procedures 
developed in this thesis resulted in successfully designing 
compensation for a system which, though not overly complex, 
was subject to somewhat severe specifications. Alternate 
methods of compensation design may have resulted in equally 
simple comepnsation but, by their trial-and-error nature, 
must be assumed to be less straightforward. It it there- 
fore concluded that the proposed compensation technique is 
valid and presents the designer with a powerful and labor- 
saving tool. 
9.2 Recommendations. Several areas for further study present 
themselves upon consideration of the results of this thesis 
and of the steps by which these results were obtained. 
Specific areas requiring further work are: 
a. More exact determination of the deviation from 
second-order response caused by the existence of 
zeros and additional roots. 
be Careful consideration of the step-by-step procedure 
to eliminate possible duplication of effort when order 
raising compensation is used. 
Areas of research which should provide beneficial results are: 


ae Examination of the means by which specifications can 
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be related to dominant third-order response. 

b. Consideration of creating non-dominant oscillatory 
modes at harmonics of the dominant mode, in order to 
better reproduce the wave form of the input. 

c. Examination of the means by which this method may 
be extended to include treatment of other than posi- 


tional servomechanisms. 
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APPENDIX A 


A DETAILED EXAMINATION OF THE EFFECT OF a ON 


OF A SECOND ORDER 


To examine the effect of a, on the deviation of the third 
order response from the second order, consider the variation 


of the two systems if a unit step is applied to each. 


Second Order: Out(s) > ern (A.1) 
s(s +2Cw Stw ) 
1 IL -ÜUo, t 2 
0ut(t) = — * —— — se "RM sinwe IK e taf) 
a Ji-C€ q e n i 
(A32) 
where 9,- Phase lag of ideal (2nd order) response - 
m? 
- arctan = 
an 
Third Order: Out(s) 9 — — —— 45 —— — —5- (4.3) 
s(sta,w,)(s +2€0, Stw, ) 


where 20, is included to make steady state values equal. 


-a0, t 
Out (t) = — - = + 
wo a [a“-2%4+1] 


a t 2 
a, € On” sinwI-X wn ttp 


a [Ol eten) 7/2 iege 
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where i= Phase lag of actual (3rd order) response - 


- arctan Lx - arctan VI s 
- ant 


Define deviation as follows: 


D = actual - ideal (A5) 
ideal 


Substitution of“equations (A.2) and (A.h) into equation 
(A.5) yields: - 


2 -(a-CO)ugt MT? 
ALS VE a TE a cp) 
[a~-2%a+1 ] [a“-2¢0+}) 


rat ¿opt + sin“ ay t e 91) 





(A.6) 


This deviation is a function of t,a, and © of the dominant 
roots. Ideally we would like to find the maximum deviation, 
no matter when it occurs. This, however, is a great under- 
taking which is beyond the scope of this paper. To simplify 
the problem we have picked the time corresponding to the peak 


overshoot of the second order response: 


vIe 
zu 





5 + arctan 
T = 
p on 


Substituting this time into equation (A.6): 





m2 
AR ae: .2 sin(Z*arctan cs = a, 
2 Uum EN m ms s 
D = [af -2ča+1] _ hem 


re ea & 1 
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m2 
but Sin (Z+arctan 1-5 ) = ak 
2 dt P MEME 


1 





so: DE AT e^ (o7. Ty + al+] (1.8) 
[14/4 -t^ e J95Tp](.9-2t a1] 

This equation is plotted in figure A.l. Figure A.2 represents 

the result of similar consideration of the effect of a zero 


at aw, on a second order response. 


= 


— 


Fxamination of these curves will yield an appropriate 
value of a to insure the desired degree of dominance of the 


second order response. 
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APPENDIX B 


THE MODEL CHARACTERISTIC EQUATIONS 


B-1 Second-Order Equations 
2 = 
S + 815 + ay = 0 


Minimum Model: 


aj 7 220, 
" 2 
ao On 
Working Model: 
a, = eto, 
E 2 
ao an 


B-2 Third-Order Equations 


3 2 e 
s + 858 + a,s + ay = 0 
Minimum Model: 
a = (2e+a, Jo, 
a. = (142064. )@ 2 
1 1 n 
" 3 
ay = (a, )o, 
Working Model: 
ay = (2¢+a, Jo, 
a, = (1+2%a.)w = 
1 ol 
de 3 
See (ao, 





90 
B-3 Fourth-Order Equations 


er + a,s" + ans * 6,9 * a, > O 
Minimum Model: 
ly eto, 
a = (1+0, Ja, * 
aj 7 (20a, *)a, 7 


Ls (0,0, * 


Working Model (Non-dominant roots complex) 


ay = (22,+20,% Ja, 

ET (Ithaky tota an 
a = (20 taerien Ja 
LN (0,0, * 


Working Model (Non-dominant roots real) 


ho $ (22+P,)0, 
a, = (1+208,+8,)0, 
a, = (B,+20B,)0, > 
"e (Pa, * 
where 
E w^ 
es 
B- Fifth-Order Equations 
> + as" T a5? + ans mae cr de = 0 


Minimum Model: 


A c DN 
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B- Fifth Order Equations (cont'd) 
Minimum Model (cont'd) 


84 7 (L+a, «2t, o, Ja," 

a, = (2€, 2,5 «,*a, 22, )a, 
A > (0,220, 0,2%, * 

A, = AY 


Working Model (One pair non-dominant roots complex) 


a, = (20,420, Cota, )w, 


w 
l 


e 2 2 
(neo e gto; er genee aa 


5 

a = (20 LMP] E Era +4a. an C. C +a 2a > 

2 IN H 3 2 l 27] 2 D OMEN 
E 2 2 o »% 

a, = (a, +20, 450,420, ant, Jo, 


2 
O (dy a5 )0, 


w 
l 


Working Model (Non-dominant roots real) 
a, = (20, +B, ao, 


2 
(14+2 8, +85 )w, 


3 
(Bj +20) 8 o+Bo wo, 
T T 
1 (P7+28,B3)0, 
5 


w 
l 


w 
tl 


(8320s 


where 
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B-5 Sixth-Order Equations: 
6 5 3 L 


g + A5" + aus" + x + "T + a,s + a, = 0 


Minimum Model: 


ag = (22,20, 
A a e e 

ap = (Itaj tas Jog 
" pr 2 38 

e A 
ME Ps 390 4 

a, = (a, "4 P" a5 Jo, 
" 2 2 

a, = (20,4, a, Ja, 

a = (a 2, m 6 

O T “ez n 


Working Model (Both pairs non-dominant roots complex) 
ag = (20, +20) 0542058, )0, 
a, = (1+ha t,t +a ? ur Co, HEC, Ca, o, FÀ RU e 
n rie" PR S Q'3 1 2*2 n 
E 2 2 
a. = (20, 05420, 0 +2050 3480, 450) 056342030) a, 


2 2 B 
420,05 +20, 05°C, )w, 


s 2 2 c 2 
a, = (oj ha oC C lia aC Coo, HC O00, 


P 211 
um eo oe 

i 2 2 2 2 B 
aj 7 (2049, «592a, 0, (,*2a4 a, (Ja, 


ae 2 a 20 6 
Y 22 n 


Working Model (One pair non-dominant roots real) 
ag = (2t *2e, 0*9. Qa. 


2 2 
a en P. rho 


fo 
ipe 
i 


E 2 2 
ay MO uina oy ec vem c TEE 


| 3 
+0,B2L 230, 
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B-5 Sixth-Order Equations (cont'd): 


Working Model (One pair non-dominant roots real) 


fb 
ti 


2 2 2 4 
(a1 *2a, B, C5*20, a4 B. Bra Boč ota Pa. 
á 2 2 5 
aj = (oe, B,+20,B,0,+20, "BA, 20, 


2 6 
m = Ba, 


f 
l 


P = faa 


Working Model (All non-dominant roots real) 


ar (20, +8, ow, 
3h 
È 3 


2 
(1420, 8, +8, a, 


w 
I 


w 
l 


h 4 
> (p5*20, B,*P, a. 


m 
l 


2 2 


w 
l 


= (Ba, 


Pa = Sy 
Ba = A a RA 
Pz = A AO o 


BL = 4,0,0,.a 


f 29 9h 


B-6 Seventh-Order Equations: 
7 6 5 


28 





Minimum Model: 


a 


"^ 


6 


=" 


2 2 
= (a, +a, +0, 0, “420, 0 a4*20, 0, 230, 
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(20,403 )0, 
(Ira, ta,“ #2014, )u, 2 


3 
(2t, «Arat es A^ 1332 53), 
a ET 


2 2 2 2 2 2 
(205%) a, ta, atan agta a, a3)0 

24 La 2 6 
unas N 1 %2 53945 


24 2 7 
(a a, a, a3 )0, 


Working Model (Some roots complex) 
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FÉ 


"i 


where Án refers to 


sixth order eqns. 


Working Model (All 


a 


ag 
au 


a 


6 


3 


(A,ta)a, 

(A,+94 50, 
(A tal, Ja, 
(AptaA a, 
(A +eA,)o, 
(A, tA, wo, 
(Ad, 


coefficients of same type model in 


non-dominant roots real) 
(26,+8, Jw, 
(1+2%,P,+B,)w, 
(B,+2%,B2+P3)0,” 


(8552€, 8*5 Da, 


5 


ls 





where 


99 


: 5 


a, = 
6 
a, = (By +20,85)0, 
, 7 
Bo E (Be Jo, 


> a. ta ta ta 


PU ño AA 


= MAD OU MEUSE EUR, e) pA yh tg Fat gM hore osa dao 


. 


e 


=z a a _a_ta_a_a, ta_a_a_ta_a_a +44 a,_a +a.a a 
234 3*5 


EO. 199» Fes 2.5 5 


DOCERE 
pec gw WOO or ASA 
2 95949, 5 
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